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基于 MODIS TVDI 和 模糊 数学 方法 的 藏 北 地 区 
旱情 等 级 遥感 监测 
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摘 要 : 为 了 实现 对 藏 北 区 域 范 上 


内 春 夏 旱情 的 动态 连续 监测 ,基于 温度 植被 干旱 指数 (TVDI) 和 模糊 数学 方法 建 


立 了 瑰 感 干旱 的 划分 标准 ,研究 时 有 段 为 1980 一 2017 年 。 首 先 利 用 MODIS 产品 数据 计算 TVDI, 然 后 根据 气象 干旱 
等 级 监测 结果 ,采用 模糊 数学 法 建立 基于 MODIS TVDI 的 干旱 等 级 划分 标准 ,并 对 监测 结果 进行 精度 验证 ,最 后 分 
析 了 近年 来 藏 北 地 区 旱情 的 时 空 变化 特征 。 得 到 的 主要 结论 :四 基于 归 一 化 植被 指数 (NDVI) 和 增强 植被 指数 


(EVI) 计算 得 到 的 温度 植被 干旱 指数 TVD 和 TVD, , 均 与 20 cm 实测 土壤 水 分 含量 在 0. 05 的 水 平 达 到 显著 相 
关 ,TVDIE 的 决定 系数 更 高 ;@) EF TVD 将 旱情 划分 为 无 旱 、 轻 旱 、 中 旱 、 重 旱 \ 特 旱 5 个 等 级 ,其 中 , 据 此 标准 获 
得 的 藏 北 地 区 旱情 等 级 与 气象 干旱 等 级 监测 结果 大 体 一 致 ;@ 近年 来 藏 北 地 区 旱情 整体 不 太 严 重 , 且 总 体 趋 缓 ,其 
中 ,2009 年 最 严重 ,发 生 中 旱 及 以 上 旱情 的 区 域 面积 达 24% ,年 内 旱情 在 6 月 最 严重 。 就 旱情 的 空间 分 布 特征 而 


ll 


,研究 区 西南 部 和 中 部 干旱 比较 严重 ,北部 和 东南 部 相对 较 轻 。 研 究 成 果 可 为 藏 北 地 区 干旱 监测 提供 数据 支撑 ， 


遥感 干旱 等 级 的 划分 方法 可 为 其 他 地 区 的 干旱 研究 提供 参考 。 
关键 词 : 温度 植被 干旱 指数 ; 模糊 数学 法 ; 干旱 等 级 ; MODIS; 藏 北 地 区 


旱灾 是 我 国 最 严重 的 自然 灾害 之 一 ,长 时 间 较 
为 严重 的 干旱 不 仅 直接 影响 农业 生产 、 社 会 经 济 ,而 
且 还 会 出 现 河 水 断 流 .绿洲 缩减 .土壤 荒漠 化 VDE 
暴 增多 等 衍生 灾害 ""。 早 灾 是 西藏 最 严重 的 农业 
气象 灾害 之 一 ” ,每 年 均 有 不 同 程度 的 发 生 , 对 农 
牧 业 生产 的 影响 极 大 ” 。 然 而 ,由 于 西藏 高 原 基础 
台 站 观测 资料 不 足 ,给 干旱 监测 带 来 许多 困难 。 在 
全 球 气候 变 暖 的 背景 下 ,极端 干旱 的 面积 不 断 扩大 ， 
西藏 西南 部 干旱 发 生 频 率 较 高 ” 。 因 此 ,开展 藏 北 
地 区 的 干旱 动态 监测 具有 重要 的 现实 意义 。 

遥感 技术 为 大 范围 .长 时 间 序 列 的 干旱 监测 提 
供 了 有 效 的 手段 5 。 大 量 研究 表明 ,温度 植被 干旱 
虽 数 (TVDI) 模型 能 较 好 地 反映 大 范围 的 土壤 水 分 
状况 。 这 种 方法 在 我 国 的 适用 性 也 得 到 了 广泛 
的 证 明 ,如 薛 天 翼 等 ”研究 表明 ,土壤 相对 湿度 和 
降水 量 数据 与 TVDI 都 呈 明 显 的 负 相关 。 近 年 来 ， 
不 少 学 者 对 TVDI 监测 土壤 水 分 进行 了 发 展 和 改 
进 ,改进 后 的 TVDI 监测 干旱 效果 更 好 。 王 行 汉 
等 基于 增强 植被 指数 对 华南 农业 旱情 进行 遥感 
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监测 ,结果 能 够 准确 地 体现 出 不 同等 级 旱情 的 空间 
区 域 分 布 状况 。 杨 秀 海 等 ” 曾 成 功 地 采用 温度 植 
被 干旱 指数 (TVDI) 评 价 了 青藏 高 原 的 旱情 , 宋 春 桥 
等 中 亦 采用 TVDI 对 于 藏 北 的 土壤 湿度 分 布 进行 了 
研究 。 

为 准确 把 握 干 旱 对 农 牧 业 生产 的 实际 影响 , 干 
旱 监测 需要 对 旱情 的 严重 程度 进行 准确 的 评估 。 基 
于 遥感 干旱 指数 TVDI 进行 旱情 等 级 的 估算 ,TVDI 
分 级 标准 人 为 性 很 大 ,从 当前 发 表 的 文献 中 可 以 看 
到 不 同 研究 中 TVDI 的 分 级 标准 不 尽 相同 呈 ) 。 目 前 
多 采用 经 验方 法 ,如 自然 断 点 法 所] 。 显 然 , 这 种 估 
算 方 法 未 考虑 到 干旱 等 级 在 区 域 中 的 适用 性 ,因此 ， 
有 必要 科学 地 建立 干旱 等 级 的 区 域 标准 。 关 于 干旱 
等 级 评估 的 标准 ,尽管 监测 方法 和 指标 众多 ,但 气象 
部 门 普遍 认为 基于 降水 距 平 获得 的 干旱 等 级 评估 结 
果 更 加 符合 实际 情况 。 因 而 ,在 气象 业务 工作 中 多 
采用 基于 降水 距 平 率 的 干旱 等 级 标准 。 为 此 ,参考 
气象 干旱 监测 的 结果 ,建立 基于 遥感 干旱 指数 的 干 
旱 等 级 划分 标准 ,并 用 于 旱情 的 监测 与 分 析 , 对 于 及 
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刘 一 哲 等 :基于 MODIS TVDI 和 模糊 数学 方法 的 藏 北 地 区 旱情 等 级 遥感 监测 


时 获取 藏 北 地 区 区 域 尺度 的 旱情 分 布 及 发 展 情况 具 
有 重要 的 意义 。 

根据 模糊 数学 的 隶属 度 理论 可 对 受到 多 种 因素 
制约 的 事物 或 对 象 作出 一 个 总 体 的 评价 。 它 具有 结 
果 清 晰 ,系统 性 强 的 特点 ,能 较 好 地 解决 模糊 的 、 难 
以 量化 的 问题 ,适合 各 种 非 确 定性 问题 的 解决 。 
exe vi 利用 规范 值 表示 的 湖泊 富 营养 化 指标 ( 即 
得 标 总 体 ) 用 同一 个 “等 效 " 规 范 指标 代替 ,继而 采 
用 模糊 数学 地 方法 成 功 地 划分 了 湖泊 富 营养 化 的 评 
价 标准 ,并 在 福州 市 山 仔 水 库 水 质 监 测 上 取得 了 很 
好 的 效果 。 在 遥感 干旱 等 级 的 评价 中 ,TVDI 可 被 视 
为 一 个 综合 的 干旱 监测 指标 ,依据 降水 距 平 获得 的 
干旱 等 级 结果 可 作为 参考 的 标准 ,利用 TVDI 建立 
干旱 等 级 的 划分 特征 满足 模糊 数学 使 用 标准 。 因 
此 ,本 文 拟 采用 模糊 数学 方法 建立 藏 北 地 区 基于 
TVDI 的 干旱 等 级 划分 标准 ,以 期 在 藏 北 地 区 建立 更 
加 符合 客观 事实 的 干旱 等 级 的 监测 方法 和 标准 ,为 
气象 部 门 更 好 地 服务 于 农 牧 业 生 产 提供 科学 依据 ， 
也 为 其 他 区 域 建立 基于 遥感 指数 的 旱情 等 级 监测 提 
供 方法 的 借鉴 。 


1 数据 基础 与 方法 


1.1 研究 区 概况 

藏 北 地 区 是 指 我 国 西藏 范围 内 介 于 78°%23'40” ~ 
95°10'46"E ,29°40'40" ~ 35°42'55"N 之 间 的 区 域 (图 
1) 。 该 区 平均 海拔 4 500 m, 气 候 寒 冷 干燥 ,年 降水 
量 247.3 ~513.6 mm。 降 水 东南 多 ,西北 少 , 且 季 节 
性 强 ,时 空 分 布 不 均 。 受 地 形 和 气候 分 异 的 影响 , 植 
被 类 型 由 东南 到 西北 依次 是 亚 高 山 玖 林 一 灌 从 草 
向 ,高寒 草 向 高寒 草 原 ,高寒 荒漠 草原 4 大 基带 及 


相 邻 的 过 渡 亚 带 。 区 内 拥有 广阔 的 天 然 草地 资源 ， 
各 类 天 然 草地 面积 占 到 全 区 总 面积 的 80% , 占 西藏 
天 然 草 地 面积 的 59% ,因而 是 西藏 重要 的 畜牧 业 生 
产 基地 。 根据 宋 春 桥 等 的 研究 成 果 , 受 水 平 
和 垂直 地 带 性 气候 的 影响 , 藏 北 地 区 水 热 和 气候 条 
件 差 异 很 大 ,因而 ,全 区 可 分 为 昆仑 高 山高 原 亚 寒带 
干旱 区 . 芜 塘 高 原 湖 盆 亚 寒带 半 干 旱 区 AA a 
寒带 干旱 区 和 那曲 高 山谷 地 高 原 亚 寒带 半 湿 润 区 4 
个 气候 区 (图 la) 。 该 区 域 的 旱情 主要 发 生 在 春 、 夏 
两 季 "“ 。 特 别 是 2009 年 人 夏 以 来 ,西藏 降水 量 偏 
少 , 气 温 持 续 偏 高 ,西藏 出 现 大 面积 中 度 旱情 ,部 分 
地 方 出 现 了 30 年 一 遇 的 气象 干旱 。 受 高 原 地 形 的 
限制 ,全 区 的 气象 站 数量 稀少 且 分 布 不 合理 (图 
1b) ,难以 为 农 牧 业 提供 准确 全 面 的 气象 干旱 数据 
SCH. 
1.2 数据 与 处 理 

使 用 MODIS 产品 组 按照 统一 算法 研发 的 
MOD13A2 全 球 1 km 分 辨 率 植 被 指数 16 d 合成 产 
品 和 MOD11A2 全 球 1 km 分 辨 率 地 表 温 度 8 d 合成 
L3 产品 。 数 据 的 时 间 为 2002—2018 年 每 年 的 3 一 8 
月 ,所 用 区 域 标识 为 h24v05、h25v05、h25v06、 
h26v05 。 其 中 ,植被 指数 产品 包含 NDVI EVI 产品 ， 
通过 NASA 提供 的 MODIS Reporjection Tool (MRT) 
读 取 文件 .拼接 .投影 ,并 通过 S-G 滤波 5 降 噪 处理 
后 生成 TIF 格式 的 植被 指数 影像 。 温 度 产品 取 平 均 
方法 合成 16 d 数据 。 

为 验证 MODIS TVDI 对 旱情 的 监测 能 力 , 分 别 
于 2015 年 8 月 和 2016 年 6 月 利用 Field Scout TDR 
土壤 水 分 计 在 实地 采集 获取 藏 北 地 区 20 cm 土壤 水 
分 含量 ,采样 点 分 布 如 图 La 所 示 。 实 地 采集 土壤 水 
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图 1 藏 北 地 区 基本 情况 


Fig.1 Basic situation in north Tibet 


于 g 


分 数据 时 ,以 1 km x1 km 为 样 方 ,每 个 样 方 按照 S 
形 间隔 , 取 10 个 土壤 水 分 实测 值 , 取 平 均值 作为 该 
采样 点 的 土壤 水 分 体积 含量 。 

藏 北 地 区 国家 站 和 自动 站 (图 1b) 的 月 降水 量 
数据 由 西藏 高 原 气 象 研究 所 提供 ,数据 格式 为 . txt 
的 文本 格式 ,时 间 为 1980 一 2015 年 ,用 于 计算 旱情 
等 级 的 参考 标准 和 评估 结果 的 验证 。 其 中 ,由 于 9 
个 国家 站 建站 时 间 较 长 , 且 为 人 工 监测 站 ,数据 准 
确 , 因 此 ,利用 国家 站 的 数据 计算 出 36 a 内 3 一 8 月 
的 月 平均 降水 量 , 再 计算 研究 时 段 内 月 降水 距 平 百 
分 率 , 据 此 获得 站 点 所 在 区 域 的 气象 干旱 等 级 。 自 
动 站 的 数据 时 间 较 得 ,剔除 异常 值 处 理 后 ,用 于 对 朋 
感 旱情 等 级 监测 的 结果 准确 性 验证 。 

1.3 研究 方法 

1.3.1 温度 植被 指数 的 计算 ” 先 用 MODIS 植被 指 
数 (NDVI) 和 增强 植被 指数 (EVI) 分 别 与 MODIS 陆 
表 温度 (LST) 建 立 TVDI, 然 后 通过 实测 20 cm 土壤 
水 分 验证 两 种 TVDI 监测 的 精度 ,选取 更 优 数据 。 
考虑 到 藏 北 地 区 植被 稀疏 .返青 期 和 生长 期 较 短 ,为 
此 ,计算 时 将 植被 指数 在 0 以 上 的 MODIS 数据 全 部 
参与 干 湿 边 拟 合 。 应 用 经 过 S-G 滤波 平滑 降 噪 处 理 
后 的 NDVI EVI 以 及 LST 预 处 理 后 的 产品 数据 , 按 
照 每 个 时 期 的 干 湿 边 方程 系数 , 反 演 求 得 所 有 像 元 
的 TVDI 值 。 

1.3.2 基于 TVDI 的 干旱 标准 划分 方法 ”考虑 到 
研究 区 内 气候 特征 的 空间 差异 和 气象 站 点 的 分 布 情 
况 , 首 先 将 研究 区 分 为 半 湿 润 区 和 半 干 旱 区 两 部 
分 5。 其 中 ,图 1 中 的 那曲 高 山谷 地 高 原 亚 寒带 半 
湿润 区 归 为 半 湿 润 区 ,其 余 3 个 气候 区 为 半 干 旱 区 。 
然后 ,采用 模糊 数学 法 对 两 部 分 分 别 建立 各 自 的 干 
旱 等 级 标准 ,基本 步骤 包括 :第 一 ,确定 被 评价 对 象 
的 因素 集 ,其 中 ,因素 集 为 刻画 评价 对 象 的 各 种 评价 
因素 ,在 这 里 是 9 个 国家 站 所 在 栅 格 的 降水 距 平 百 
分 率 等 级 和 TVDI 值 。 为 了 使 干旱 标准 具有 代表 
PE, AS SCHR 2002—2014 年 的 降水 距 平 百分率 和 TV- 
DI 平均 值 用 于 构建 因素 集 ; 第 二 ,确定 评语 集 ,评语 
集 是 评价 者 对 被 评价 对 象 可 能 作出 的 各 种 总 的 评价 
结果 组 成 的 评语 等 级 的 集合 ,依据 气象 干旱 等 级 的 
标准 ,评语 集 为 无 旱 、 轻 旱 、 中 旱 、 重 旱 、 特 旱 ; 第 三 ， 
将 TVDI 干旱 分 级 特征 向 量 归 一 化 确定 评价 对 象 的 
权重 向 量 : 


(1) 


式 中 :Xi 表示 归 一 化 后 第 i 个 特征 向 量 值 ;X; 表示 
第 i 个 特征 向 量 值 ,i=1,2,…,n。 本 人 研究 中 n=4。 

其 次 ,依据 不 同 干旱 指标 等 级 的 常规 分 级 标准 ， 
进行 各 类 指标 评价 函数 设计 ;然后 ,根据 评价 函数 以 
及 评价 矩阵 ,计算 模糊 隶属 度 , 即 得 到 模糊 关系 矩阵 
Ro), 


Phyo Tig Ty 
ra Ta 72 

R= (2) 
Tal Tn wl Tan 


式 中 ;7, 表 示 某 个 被 评价 对 象 从 综合 指标 un (TV- 
DI) RAM SRB LA vw 的 隶属 度 , 这 里 vw WA 
象 干旱 等 级 。 

最 后 ,将 模糊 关系 矩阵 与 因素 的 权重 矢量 进行 
模糊 运算 ,得 到 模糊 综合 评价 结果 。 

B=XxR (3) 

式 中 :B 为 得 到 的 TVD 等 级 值 ;X 为 归 一 化 确定 评 
价 对 象 的 权重 向 量 集 。 
1.3.3 干旱 等 级 划分 标准 的 检验 利用 研究 区 23 
个 自动 站 (图 1b)2010 一 2017 年 的 数据 计算 气象 降 
水 距 平 率 评估 干旱 等 级 ,并 与 其 所 在 栅 格 的 同期 遥 
感 干 旱 等 级 监测 结果 进行 对 比 ,以 验证 遥感 旱情 等 
级 监测 结果 的 正确 性 。 需 要 指出 的 是 ,由 于 研究 区 
北部 的 昆仑 高 山高 原 亚 寒带 半 干 旱 区 和 东北 部 的 青 
南 高 原 亚 寒带 半 干 旱 区 缺少 气象 站 点 数据 ,其 监测 
结果 不 能 准确 订正 。 但 由 于 两 个 区 域 的 气候 特征 与 
芜 塘 高 原 湖 盆 亚 寒带 半 干 旱 区 中 北部 和 东北 部 相 
似 , 地 理 位 置 上 前 者 距 狮 果 河 \ 改 则 气象 站 点 较 近 ， 
后 者 与 双 湖 、 扎 仁 气象 站 相近 , 故 认为 本 研究 所 建立 
的 半 干 旱 区 的 干旱 等 级 划分 标准 适用 于 该 区 域 。 另 
外 ,那曲 高 山谷 地 高 原 亚 寒带 半 湿 润 区 以 及 羌 塘 高 
原 湖 倪 亚 寒 带 半 干 旱 区 南部 的 普兰 县 所 在 区 域 ,由 
于 降水 较 多 ,气候 相对 湿润 , 故 采用 半 湿 润 区 的 干旱 
等 级 划分 标准 。 


结果 与 分 析 


2.1 TVDI 模型 的 建立 与 验证 

2.1.1 植被 指数 一 温度 特征 空间 图 2 为 2017 年 
3 一 8 月 每 个 月 第 2 期 NDVI-LST 与 EVI-LST 特征 空 
间 与 干 湿 边 方程 ,由 图 2 可 以 看 出 ,研究 区 植被 指 
数 一 地 表 温 度 特 征 空间 的 特点 :四 特征 空间 呈现 出 


2 
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注 :7; 为 像 元 地 表 温 度 值 。 


三 角形 , 干 边 和 斜率 小 于 0, 说明 随 着 植被 的 增加 ,最 
大 值 呈 减 小 趋势 ;@ 当 NDVI 数值 超过 某 一 


图 2 2017 年 MODIS NDVI-LST 和 EVI-LST 特征 空间 
Fig.2 Feature space of MODIS NDVI-LST (left) and EVI-LST (right) in 2017 


时 , 近 地 表 温度 的 最 大 值 随 着 NDVI 数值 的 增加 不 
BE 。 断 降低 ,而 最 小 值 却 不 断 增加 ,因此 , 干 边 方程 拟 合 


干 蛙 t 


效果 优 于 湿 边 方程 ;8) NDVI 和 EVI 两 种 数据 源 的 
散 点 分 布 大 体 相同 ,但 NDVI 范围 主要 集中 在 0 ~ 
0.8,EVI 范围 则 集中 在 0 ~0.6;@ 基于 增强 植被 指 
数 (EVI) 得 出 的 拟 合 系数 高 于 同时 期 植被 指数 ND- 
VI, AMEE EVI 数据 拟 合 结果 精度 较 好 。@) 在 
3 月 4 月 和 5 月 植物 生长 初期 ,LST 与 EVI 存在 正 
相关 关系 ,6 月 7 月 和 8 月 时 ,LST 与 EVI 相关 关系 
为 负 。 


平 率 等 级 TVDIs| ;这 里 的 降水 距 平 率 等 级 是 根据 
国家 站 6 一 8 月 的 月 平均 降水 量 数据 计算 得 到 ,等 级 
标准 根据 《GBT 20481 一 2006 气象 干旱 等 级 》 划 分 为 
无 旱 、 轻 旱 、 中 旱 、 重 旱 和 特 旱 。TVDI: 为 气象 数据 
同期 站 点 所 在 栅 格 的 TVDI: 平均 值 。 为 了 使 模型 
具有 通用 性 ,对 2002—2014 年 每 年 得 到 TVDI 求 平 
均值 作为 最 终 评价 和 矩阵 的 结果 。 

(2) 确定 评价 对 象 评语 集 V= | 无 旱 、 轻 旱 、 


中 


在 植物 生长 受 能 量 或 温度 作为 主要 限制 因子 

时 ,LST 与 EVI 呈正 相关 关系 ,在 植物 生长 受 水 分 作 
为 主要 限制 因子 时 ,LST 与 NDVI 呈 负 相关 关系 2 。 
故 本 文选 取 6 一 8 月 对 藏 北 和 干旱 进行 研究 。 
2.1.2 TVDI 模型 的 验证 性 评价 ”根据 前 人 研究 
TVDI 与 土壤 水 分 含量 的 验证 结果 ,TVDI 与 20 cm 
土壤 水 分 的 相关 性 最 佳 ” 。 据 此 ,本 文 利用 2015 
年 和 2016 年 实测 20 cm 土壤 水 分 与 同期 反 演 得 到 
的 TVDI 进行 相关 性 分 析 ,以 检验 TVDI 对 土壤 含水 
量 的 表达 能 力 。 

分 别 以 TVDI、( 基 于 NDVI 计算 得 到 ) 和 TVDT'， 
(基于 EVI 计算 得 到 ) 为 横 坐 标 ,实测 20 em 土壤 水 
分 为 纵 坐 标 , 制 成 散 点 图 (图 3)。 从 图 3 中 可 以 看 
出 ,两 种 指数 反 演 的 TVDI 与 实测 20 cm 土壤 水 分 
的 相关 性 都 达到 了 0. 05 的 显著 性 标准 。 其 中 ， 
TVDI, 与 实测 20 cm 土壤 水 分 的 相关 系数 高 于 TV- 
Diy 与 实测 20 cm 土壤 水 分 的 相关 系数 。 该 结果 进 
一 步 反 映 出 基于 增强 植被 指数 (EVI) 与 地 表 温 度 计 
算得 到 的 TVDT 更 适用 于 藏 北 地 区 的 干旱 监测 。 
2.2 基于 模糊 数学 方法 的 干旱 等 级 划分 
2.2.1 干旱 等 级 划分 标准 的 建立 ”采用 模糊 数学 
的 方法 ,建立 基于 TVDI; 干旱 等 级 的 标准 ,具体 过 
FEU F: 

(1) 确定 评价 对 象 因素 (指标 ) 集 U = | 降水 距 


(a) TVDI, 


y=-34.57x+38.2 


wor a 
© © 


土壤 水 分 /% 
= 


0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
TVDI 


F ER RH. 

(3) 利用 降水 距 平 百分率 按照 《GBT 20481— 
2006 气象 干旱 等 级 》 得 到 5 级 干旱 等 级 的 气象 干旱 
分 级 结果 ,再 统计 同一 分 级 结果 对 应 的 所 有 TVD, 
均值 。 考 虑 到 热量 是 高 海拔 地 区 春季 植物 的 生长 限 
制 因子 ,此 时 LST 与 EVI 存在 正 相 关 关 系 ，TVDI 不 
适合 干旱 监测 , 故 只 选取 夏季 6 一 8 月 每 个 月 按 
降水 等 级 统计 对 应 TVDI, 均值 ,最 后 组 成 评价 矩阵 
( 表 1)。 

(4) 根据 模糊 数学 的 方法 , 设 和 = (x ,x,,…， 
Xm ) 为 权重 ( 权 数 ) 分 配 模 糊 矢 量 , 其 中 表示 第 i 
个 因素 的 权重 ,要 求 x >0,Yx =1。 这 里 的 因素 权 
重 指 常规 TVDI 干旱 分 级 标准 0. 40 .0. 60 .0. 75 、 
0.86) ,将 它们 归 一 化 后 得 到 评价 对 象 的 权重 向 量 
X = (0.15 .0.23 .0.29 .0.33) 。 

(5) 指标 评价 函数 的 设计 如 表 2 所 示 , 其 中 , 建 
立 2 个 气候 区 的 标准 时 采用 相同 的 评价 函数 。 

(6) 根据 评价 函数 构建 评价 和 矩阵 ,计算 模糊 隶 
属 度 , 即 得 到 模糊 关系 矩阵 ,分 别 获得 半 湿 润 区 和 半 
干旱 区 2 个 模糊 隶属 度 和 矩阵 ( 表 3) 。 

(7) 计算 得 到 阔 值 向 量 ,建立 藏 北 地 区 MODIS 
TVDL 干旱 等 级 划分 标准 ,具体 结果 如 表 4 所 示 。 
2.2.2 标准 的 适用 性 讨论 与 验证 气象 干旱 是 指 
降水 与 蒸 散 收 支 不 平衡 造成 的 异常 水 分 短缺 现象 ， 


y=-73.99x+66.2 
: R°=0.607 9 


土壤 水 分 /% 
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TVDI 


图 3 TVDIy 和 TVDI: 与 实测 20 cm 土壤 水 分 散 点 图 
Fig.3 TVDIN and TVDI, with a 20 - cm-deep soil water dispersion plot 
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表 1 评价 矩阵 
Tab.1 The evaluation matrix 
ai 半 湿 润 区 六 干旱 区 
无 旱 轻 旱 中 旱 重 早 无 时 干旱 中 旱 重 旱 
6 0.25 0.42 0.60 0.74 0.58 0.65 0.78 0.88 
7 0.38 0.58 0.66 0.79 0.61 0.69 0.82 0.91 
8 0.41 0.40 0.64 0.75 0.64 0.60 0.74 0.85 
R2 各 指标 评价 函数 表 
Tab.2 Evaluation functions of the indexes 
评价 指标 评价 函数 取 值 范围 评价 指标 评价 函数 取 值 范围 
时 0.40 ils 0 <TVDI, <0. 40 时 9:13 is 0 <TVDI,; <0.75 
ae TVDI, es Pe TVDI, -0.60 BARSA 
TVDI,, -0.40 TVDI -0.75 
— 一 -一 TVDI,, >0. 40 —_* TVDI >0. 75 
TVDI, TVDI,, -0.60 
时 :00 OI 0 <TVDI, <0.60 重 早 peas 0 <TVDI, <0. 86 
Te TVDI -0.40 pan a TVDI,, -0.75 _ 
TVDI,, -0.60 TVDI, -0.86 
tt E TVDI, >0. 60 eee ate TVDI, >0. 86 
TVDI,, -0.40 TVDI,, -0.75 
表 3 模糊 隶属 度 和 矩阵 
Tab.3 Fuzzy membership matrix 
AOD 半 湿 润 区 半 干 旱 区 
无 旱 轻 旱 中 旱 重 早 无 早 旱 中 旱 重 旱 
1 0.28 0.47 0.64 0.68 0. 43 0.57 0.76 0. 82 
2 0.34 0. 42 0.56 0.62 0. 62 0.69 0.80 0. 87 
3 0.43 0. 66 0.79 0.81 0. 64 0.74 0.85 0.92 
#4 + MODIS TVDI, 的 干旱 等 级 监测 划分 标准 
Tab.4 Classification standard of drought based on the MODIS TVDI: 
等 级 ”类 型 半 湿 润 区 EP RR 等 级 AA 半 湿 润 区 半 干 旱 区 
1 无 旱 TVDI <0. 37 TVDI, <0. 52 4 me 0.69 <TVDT <0. 83 0.82 <TVDI <0. 91 
2 轻 旱 0.37 < TVDI; <0.55 0.52 <TVDI <0. 71 5 特 旱 TVDI 20. 83 TVDI 20. 91 
3 中 时 0.55 < TVDI; <0. 69 0.71 < TVDI; <0. 82 


主要 研究 天 气 的 干 湿 程 度 。 农 业 干 旱 是 指 作 物 生长 
过 程 中 因 供 水 不 足 ,阻碍 作物 正常 生长 而 发 生 水 量 
供需 不 平衡 的 现象 ~“ 。 本 文 所 用 的 TVDI 数据 不 是 


TVDI, 进行 平均 值 合成 后 ,以 23 个 自动 站 为 中 心 周 
围 9 个 栅 格 值 取 平均 后 ,根据 表 4 标准 计算 遥感 干 
旱 等 级 ,再 将 结果 与 依据 降水 距 平 率 计 算得 到 的 干 


单 时 次 数据 ,而 是 半月 合成 数据 。 根 据 薄 燕 飞 5 研 
究 的 不 同 土地 类 型 下 TVDI 与 气象 干旱 指数 的 关系 


早 等 级 (气象 干旱 等 级 ) 进行 比较 。 图 4 为 6 一 8 月 
2 种 监测 结果 等 级 差 值 。 从 图 4 可 以 看 出 :遥感 干 


中 ,TVDI 与 10 cm 土壤 水 分 .TVDI 与 20 cm 土壤 水 
分 .TVDI 与 SPI( 标 准 化 降水 指数 ) ,TVDI 与 MI( 相 
对 湿润 度 指数 ) 分 别 在 对 应 多 上 相关 程度 最 高 。 同 
时 ,干旱 常用 土壤 含水 量 作为 判定 指标 ,结合 实测 


时 等 级 与 气象 干旱 等 级 多 数 情 况 下 是 一 致 的 。 夏 季 
(6 一 8 月 ) 遥 感 干旱 等 级 与 气象 干旱 等 级 监测 结果 
一 致 的 比例 均 为 85.3% ~89.7% 。 每 月 2 种 评价 
等 级 相差 1 级 的 比例 均 大 于 95% 。 由 于 当 著 歼 措 


20 cm 土壤 水 分 与 同期 TVDI 相关 性 分 析 ( 图 3) 可 
知 ,采用 TVDI 指数 在 藏 北 地 区 进行 旱情 的 监测 是 
可 行 的 。 

在 ArcGIS 中 将 2015 一 2017 年 每 月 2 期 的 


气象 站 地 理 位 置 距 离 湖 面具 有 20 m, 对 应 的 卫星 数 
据 中 一 定 包 含 了 湖面 的 信息 ,不 适合 与 卫星 数据 对 
比 , 故 做 剔除 处 理 。 以 上 结果 表明 , 基于 MODIS 
TVDL 及 干旱 等 级 划分 标准 的 遥感 干旱 监测 方法 适 
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用 于 藏 北 地 区 区 域 范围 内 的 干旱 等 级 监测 ,其 监测 
结果 与 气象 干旱 监测 结果 基本 相符 。 
2.3 旱情 遥感 监测 结果 的 时 空 动态 特征 
2.3.1 旱情 的 时 序 特征 

(1) 不 同等 级 旱情 的 年 际 变 化 特征 

统计 2002—2017 年 各 年 份 6 一 8 月 发 生 中 旱 及 
以 上 等 级 旱情 的 面积 比例 ,制作 中 旱 及 以 上 年 际 变 
化 比例 柱状 图 (图 Sa) ,可 以 看 出 :QD 藏 北 地 区 发 生 
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早 灾 的 情况 比较 常见 ,平均 每 年 有 20% 左右 的 地 方 
会 发 生 中 旱 及 以 上 旱情;@) 2009 年 灾情 最 为 严重 ， 
发 生 中 旱 及 以 上 灾情 的 区 域 面积 达到 24% ,根据 西 
藏 自治 区 防汛 抗旱 指挥 部 的 监测 结果 ,2009 年 西藏 


全 区 遭遇 了 10 a 不遇 的 干旱 ,这 与 本 研究 的 监测 评 
价 结果 一 致 。 其 次 ,2002 年 和 2015 年 的 旱情 也 相 
对 较 重 ;@ 2016 年 和 2017 年 研究 区 发 生 中 旱 以 上 
旱情 的 面积 比例 相对 较 低 ,为 11% ~12% 。 其 次 ， 
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4 6 一 8 月 自动 站 遥感 干旱 监测 与 气象 干旱 等 级 差 什 


Fig.4 Difference of drought level between the remote sensing results monitored by the automatic station and the meteorological 


drought from June to August 
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图 5 中 旱 及 以 上 面积 比例 与 各 等 级 干旱 面积 比例 


Fig.5 Proportions of moderate and above drought area and the drought area at different levels 


2011 年 的 比例 也 较 低 , 为 12.72% 。 

从 各 类 旱情 发 生 面 积 比 例 的 多 年 变化 情况 (图 
5b) 可 以 看 出 , 藏 北 地 区 的 旱情 以 轻 旱 和 中 旱 居 多 ， 
重 旱 和 特 旱 发 生 的 面积 比例 相对 较 小 。 在 严重 干旱 
的 年 份 , 重 旱 和 特 旱 的 面积 比例 明显 升 高 。 

(2) 不 同等 级 旱情 的 年 内 特征 

统计 2002 一 2017 年 每 年 6 一 8 月 同期 不 同 干旱 
等 级 遥感 监测 面积 比例 的 平均 值 ,制作 成 柱状 图 ,从 
图 6 可 以 看 出 :研究 区 在 6 月 发 生 较 严 重 的 旱情 ,发 
生 中 旱 及 以 上 旱情 的 面积 均 在 20% 以 上 ,对 应 为 夏 
量 。 夏 早 持续 到 7 月 , 随 着 降雨 增多 开始 逐渐 缓解 。 
分 析 6 月 干旱 严重 的 主要 原因 ,可 能 在 于 夏季 高 原 
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图 6 不 同 旱情 等 级 的 面积 比例 构成 及 变化 


Fig.6 Composition and change of the area proportions 


with different drought levels 
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温度 升 高 ,蒸发 量 加 大 ,加 之 主 降 水 期 还 未 到 来 , 因 
而 相对 干旱 。7 月 和 8 月 受 高 原 低 压 的 夏季 降水 影 
响 , 中 旱 及 以 上 的 比例 分 别 降 到 16. 32% 和 
12.51% ,整个 区 域 的 旱情 得 到 缓解 。 
2.3.2 旱情 的 空间 分 布 规律 ”将 2002 一 2017 年 每 
个 像 元 发 生 中 旱 及 以 上 的 次 数 比例 统计 作 图 (图 
7) ,可 以 更 加 直观 地 分 析 研 究 区 旱情 发 生 的 空间 分 
布 规律 。 可 以 看 出 , 常 受 较 严重 干旱 影响 的 地 区 主 
要 集中 在 研究 区 的 西南 部 和 中 部 ,研究 区 北部 和 东 
南部 受 旱 灾 影 响 相 对 较 轻 ,这 与 藏 北 东部 地 区 降水 
较 多 有 关 。 应 该 指出 的 是 ,由 于 半 湿 润 区 的 遥感 干 
旱 划 分 标准 低 于 半 干 量 区 ,所 以 在 半 湿 润 区 西 侧 旱 
情 略 重 于 东 侧 。 

以 旱情 最 严重 的 2009 年 为 例 ,制作 连续 的 干旱 


等 级 监测 结果 图 以 直观 地 显示 研究 区 旱情 在 空间 上 
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的 发 展 过 程 (图 8)。 从 图 8 可 以 看 出 ,研究 区 旱情 
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注 :绿色 代表 16 a 内 发 生 中 旱 及 以 上 旱情 所 占 比 例 在 0 ~ 20% 之 间 ， 
红色 代表 比例 在 80% ~ 100% 。 

7 2002—2017 年 藏 北 发 生 中 旱 及 以 上 的 次 数 比例 分 布 
Fig.7 Frequencies of moderate and above droughts in north 


Tibet From 2002 to 2017 
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8 2009 年 6 一 8 月 藏 北 地 区 旱情 分 布 特征 
Fig.8 Distribution of drought in north Tibet from June to August 2009 
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的 发 生 大 体 从 中 部 和 西部 开始 ,逐渐 扩大 影响 范围 
后 ,从 东部 开始 减弱 。 旱 情 最 为 严重 的 是 6 月 2 期， 
中 旱 及 以 上 所 占 比 例 达 到 40.8% 。 根 据 中 国 气象 
局 全 国 干旱 监测 资料 显示 ,2009 年 4 月 1 日 至 6 月 
19 日 ,西藏 中 部 及 东南 部 降水 量 为 近 10 a 来 同期 最 
少 ,特别 是 进入 6 月 以 来 ,西藏 中 南部 区 域 较 常 年 俩 
D6 成 以 上 ,符合 监测 结果 。 旱 情 加 重 时 ,人 研究 区 中 
北部 随 之 发 生 一 定 程度 的 旱情 ,主要 为 中 旱 ,不 过 该 
区 域 干 旱 影响 时 间 较 短 ,旱情 能 随 降 水 增加 较 快 消 
退 。 根 据 西 藏 自治 区 气象 局 的 气象 站 点 降水 数据 ,7 
月 和 8 月 藏 北 地 区 降水 量 较 之 前 月 份 显著 增加 , 特 
别 是 中 东部 的 那曲 地 区 ,因此 7 月 以 后 旱情 开始 由 
东 向 西 减轻 ,东部 和 北部 地 区 逐渐 恢复 正常 。 


3 结论 与 讨论 


结论 
(1) 根据 MODIS 产品 数据 获得 TVDIN 和 
TVDIs 反 演 结果 ,可 以 很 好 地 适用 于 藏 北 干旱 监测 
与 20 cm 土壤 水 分 的 相关 性 分 析 证 明 ,TVDIs XF 


3.1 


早 监测 效果 更 好 。 

(2) 采用 模糊 数学 方法 ,考虑 气候 区 差异 ,在 藏 
北 地 区 建立 了 MODIS TVDI 的 干旱 等 级 划分 标准 。 
基于 参照 气象 上 降水 距 平 率 的 气象 干旱 等 级 划分 标 
准 , 遥 感 干 旱 等 级 也 划分 为 无 旱 . 轻 旱 . 中 旱 、 重 旱 、 
特 旱 5 个 等 级 。 其 中 , 半 湿 润 区 5 个 等 级 的 划分 标 
准 依次 为 (0 ~ 0.37)、(0.37 ~ 0. 55), (0. 55 ~ 
0.69) (0.69 ~0. 83) (0.83 ~1); 半 干旱 区 5 个 等 
级 的 划分 标准 依次 为 (0 ~0.52)、(0.52 ~0.71)、 
(0.71 ~0.82) (0.82 ~0.91) (0.91 ~1) 。 遥 感 干 
旱 等 级 和 气象 干旱 等 级 的 结果 对 比 表明 ,本 研究 建 
立 的 干旱 等 级 监测 方法 和 标准 适用 于 藏 北 地 区 , 根 
据 气 候 区 建立 不 同 干旱 等 级 标准 是 必要 的 ,遥感 监 
测 的 结果 与 气象 干旱 监测 结果 大 多 情况 下 一 致 ,个 
别 区 域 和 时 间 范 围 内 存在 1 个 等 级 差 , 极 少 出 现 2 
个 或 以 上 等 级 差 。 

(3) 藏 北 地 区 的 旱情 以 轻 旱 和 中 旱 为 主 , 重 旱 
和 特 旱 发 生 的 面积 比例 相对 较 少 。 根 据 中 旱 及 以 上 
干旱 发 生 的 频次 分 布 图 ,干旱 严重 地 区 主要 集中 在 
研究 区 的 西南 部 和 中 部 。 从 干旱 发 生 的 时 间 过 程 来 
看 ,年 内 6 月 旱情 最 严重 ,7 月 中 后 期 旱情 开始 组 
解 。2002 一 2017 年 ,中 旱 及 以 上 旱情 发 生 比 例 在 
2009 年 达到 最 高 ,旱情 最 严重 。 从 旱情 频 发 的 区 域 
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来 看 ,主要 旱情 的 发 生 从 中 部 、 西 部 开始 ,从 东部 结 
Ro 
3.2 讨论 

根据 西藏 自治 区 防汛 抗旱 指挥 部 的 气象 干旱 监 
Mj ,2009 年 6 月 以 来 降水 偏 少 ,遭遇 10 年 一 遇 的 旱 
灾 ,而 在 6 月 未 西藏 部 分 地 区 迎 来 降雨 ,旱情 开始 组 
解 , 这 与 本 文 的 旱情 演变 规律 一 致 。 但 基于 TVDI 
的 干旱 监测 对 干旱 发 生 的 空间 分 布 更 详细 ,很 适合 
西藏 缺少 台 站 的 地 区 旱情 监测 。 本 研究 与 柳 锦 宝 
等 "基于 TVDI 的 遥感 干旱 监测 相 比 ,使 用 EVI 代 
替 NDVI 得 到 的 拟 合 系数 更 高 ,等 级 监测 结果 更 贴 
近 气 象 部 门 的 服务 要 求 ,能 更 好 地 为 防 灾 减 灾 决 策 
提供 依据 。 
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Abstract: A dynamic and continuous monitoring of spring and summer drought in north Tibet during the period 
of 1980 - 2017 was carried out. The temperature vegetation drought index (TVDI) , fuzzy mathematics and moni- 
tored results of meteorological drought level were used to develop the classification standard of remote sensing 
drought ,and the accuracy of monitored results was verified. The spatiotemporal variation of drought in north Tibet in 
recent years was analyzed. The main conclusions are as follows:(D TVDI, and TVDI, calculated based on the nor- 
malized vegetation index (NDVI) and enhanced vegetation index (EVI) were significantly correlated with the soil 
moisture content at 20 cm in depth at 0.05 significance level. The determination coefficient of TVDI, was higher; ®© 
Based on TVDI,, drought could be divided into four levels including without drought, slight drought, moderate 
drought , heavy drought and severe drought. On which it was obtained that the drought level in north Tibet was con- 
sistent with the monitored results of meteorological drought level ;(3) Drought in the study area was not very serious 
in recent years,and the holistical drought situation was somewhat alleviative. The drought was the most serious in 
2009 ,and the area proportion of moderate and above drought was as high as 24%. Monthly, the drought was the 
most serious in June. In terms of the spatial distribution, the drought in the southwest and the central part of the 
study area was relatively serious ,and it was relatively slight in the north and southeast. The research results could be 
referred in monitoring drought in north Tibet. 
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